Laser cooling with electromagnetically induced transparency: 
Application to trapped samples of ions or neutral atoms 
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A novel method of ground state laser cooling of trapped atoms utilizes the absorption profile 
of a three (or multi-) level system which is tailored by a quantum interference. With cooling 
rates comparable to conventional sideband cooling, lower final temperatures may be achieved. The 
method was experimentally implemented to cool a single Ca^ ion to its vibrational ground state. 
Since a broad band of vibrational frequencies can be cooled simultaneously, the technique will be 
particularly useful for the cooling of larger ion strings, thereby being of great practical importance 
for initializing a quantum register based on trapped ions. We also discuss its application to different 
level schemes and for ground state cooling of neutral atoms trapped by a far detuned standing wave 
laser field. 

PACS: 



I. INTRODUCTION 



The emergence of schemes that utilize trapped ions 
or atoms for quantum information, and the interest in 
quantum statistics of ultra cold atoms, have provided re- 
newed interest for laser cooling techniques. Starting with 
the case of a single atom, the present goal is now to laser 
cool larger numbers of atoms and to prepare them for 
applications e.g. as a quantum register. The cooling of 
a large number of particles using lasers is a prerequisite 
for coherent control of atomic systems. 

Trapped ions in linear Paul traps are currently consid- 
ered to be promising candidates for a scalable implemen- 
tation of quantum computation |l|]. Long lived internal 
states serve to hold the quantum information (qubits) 
and are manipulated coherently (single-bit gates) by laser 
beams focused on the ions [|| . The excitation of the ion's 
common vibrational motion (gate mode) provides the 
coupling between qubits which is necessary for two-bit 
quantum gate between ions. The Cirac-Zoller proposal 
for a two-bit gate [g[ requires that initially the gate mode 
is optically cooled to the ground state. Later Molmer 
and Sorensen proposed a two-bit gate which releases this 
condition and only requires cooling into the Lamb-Dicke 
regime, typically with a mean vibrational quantum num- 
ber m<l 1,1 . 

In general, the motion of a string of N cold trapped 
ions in a linear Paul trap is described by 3A^ harmonic 
modes of vibration. Cooling all modes close to the vibra- 
tional ground state is a prerequisite to realizing quan- 
tum gates with high fidelity. In typical experimental 
situations this requires cooling techniques which lead to 
lower mean vibrational state than simple Doppler cool- 
ing. Sideband cooling has been shown to achieve sub- 
Doppler temperatures [pl-p| but proves increasingly inef- 
ficient for cooling many modes. The recently proposed 
method of EIT-cooling provides parallel multi-mode cool- 



ing while its experimental effort is significantly reduced. 
The first part of this paper (Sec. pB| ) describes the 
principle of the EIT-cooling method and illustrates its re- 
sults with several realistic examples, the cooling of a sin- 
gle IIg+ ion held in a Paul trap, and of neutral Rb atoms 
confined in a far-detuned optical dipole trap (sect. p^ID| ). 
In section ID we report the first experimental demonstra- 



tion of the EIT-method, achieving ground state cooling of 
a single trapped Ca+ ion. The detailed description of all 
necessary experimental ingredients and measurements is 
meant to provide a recipe for implementing EIT-cooling 
in various experimental situations. We conclude by a 
proposal for cooling a linear string of ions (Sec. |^), and 
discuss the relevance for quantum information process- 
ing. 



II. THE PRINCIPLE OF EIT-COOLING 

A. Cooling basics 

The theory of laser cooling of atoms in a harmonic trap 
has been extensively discussed in the literature, see llQ] 
for a comprehensive review. In many experimental situa- 
tions a simple, intuitive rate equation picture |11| applies 
which can also be derived from the general theory in the 
appropriate limits pOl. In this section we recall those 
intuitive considerations, a detailed theoretical derivation 
being beyond the scope of this paper. 

Consider a trapped 2-level atom with states \g) and 
|e), excited below saturation by a laser detuned from 
resonance by A. The scattering rate, i.e. steady-state 
excited state population times decay rate, is denoted 
by VF(A). The atom is coupled to one vibrational de- 
gree of freedom, a harmonic oscillator with frequency z/ 
and quantum states |n), such that the internal transi- 
tion \g) ^^ |e) splits up into transitions \g,n) ^^ |e,n'). 



The coupling between different motional states is me- 
diated by the recoil of a laser or spontaneous photon. 
We restrict ourselves to the so-called Lamb-Dicke regime 
1] = k^J {x^) ^ 1 where the spatial extension of the mo- 
tional wave packet y/{x'^) is smaller than the laser wave- 
length A = 27r/fc. In this limit the relevant transitions 
are \n) — > \n) (carrier) and \n) — > |n ± 1) (sidebands). 



carrier absorption, being followed by either cooling or 
heating emission, causes a diffusion of the energy which 
in average contributes to heating. Therefore for opti- 
mum cooling, absorption on the cooling sideband should 
be maximized while carrier and heating-sideband absorp- 
tion should be minimized. This leads to the principle of 
EIT-cooling. 



See Fig. |i|a for an illustration. 
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order in ry. The relative probabilites of the individual pro- 
cesses are indicated, (b) Illustration of the rate coefficients 
for cooling and heating, after |lOt|. 

The probability that scattering of a photon takes the 
atom from \g,n) to \g,n± 1) has two contributions, cor- 
responding to the two possible intermediate states \e,n) 
or |e,n ± 1). The rate coefficients A± for cooling and 
heating transitions (sec Fig. |jb) arc therefore given by 



A± == rf (VK(A) + W{A T J^)) 



(1) 



Summing up all processes in lowest order 77, a rate equa- 
tion for the populations Pnit) of the states \n) is obtained 
which leads to a simple dynamical equation for the mean 
vibrational excitation {n{t)) ~ '^npn{t), 



dt 



{n) ^ -{A^ - A+){n) + A+ 



(2) 



In the more general case that the laser k vcktor is at an 
angle 9 to the direction of vibration, and that the spon- 
taneous recoils are spatially distributed with an average 
projection a on the vibrational axes ||l^, Eq. (|l|) changes 
into 



A, 



V 



{aW{A) + cos^{d)W{A t i^)) 



(3) 



Using Eqs. (0,0), the cooling limit ifi = (71(00)) = 
A+/(y4_ — A+) and the cooling time constant, TcooI = 
(A_ — A+)^^ can be easily evaluated. One finds the 
well-known cooling limits fhi/ ~ 7/2 for Dopplcr cooling 
where j > v and A = —7/2, and 771 <C 1 for sideband 
cooling where 7 < ^ and A = — j/ pO| . 

It is useful to look at the transitions displayed in 
Fig. n^a from a different viewpoint: Absorption on one of 
the sidebands causes a drift of the motional energy while 



B. EIT-cooling principle 



Electromagnetically induced transparency (FIT), like 
coherent population trapping or " dark resonance" , is a 
manifestation of quantum interference between atomic 
transition amplitudes. A review can be found in p3[ . 
FIT arises in three- (or multi-) level systems and con- 
sists in the cancellation of the absorption on one tran- 
sition induced by simultaneous coherent driving of an- 
other transition. Consider a 3-level atom as shown in 
Fig. 1, with ground state \g), stable or metastable state 
\r) and excited state |e). State |e) has hnewidth 7 and 
is coupled to \g) and \r) by dipole transitions which are 
laser-excited. FIT arises when the detunings of the two 
lasers from state |e) are equal: The system evolves into a 
coherent superposition of \g) and |r), and light scattering 
seizes. 

In order to use this situation for cooling, the transi- 
tion \r) -^ \e) is excited by an intense ("coupling") laser 
(fir ~ 7) at detuning A^ above resonance. Then the ab- 
sorption spectrum on the transition \g) -^ \e) seen by the 
other (" cooling" ) laser is described by a Fano-like profile 
|1J| which has three characteristic features, see Fig. la,b: 
the broad resonance at Ag ~ 0, the dark resonance (FIT) 



at Ag = A^ 
where 



and the narrow resonance at An 



<5=(VA2 + r!2-|A,|)/2 



is the ac-Stark shift created by the coupling laser. 



(4) 
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FIG. 2. (a) Levels and transitions of the EIT-cooling 
scheme. The inset shows schematically the absorption rate on 
\g) -^ |e) when the atom is strongly excited above resonance 
on |r) — > |e). (b) Absorption of cooling laser around Ag — Ar 
(solid line); probabilities of carrier {\n) —^ \n}) and sideband 
(|n) — > |n ± 1)) transitions when Ag — Ar are marked by 
dashed lines, (c) Dressed state picture: the cooling laser ex- 
cites resonantly transitions from \g, n) to the narrow dressed 
state denoted by \ + ,n — 1) which preferentially decays into 
\g,n-l). 
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For EIT-cooling; the laser frequencies are set to the 
dark resonance condition, Ag = A,.. Then, taking into 
account the harmonic motion, all \g,n) -^ |e,n) transi- 
tions are cancelled. Furthermore, by choosing a suitable 
Rabi frequency Qr, the spectrum is designed such that 
the \g, n) -^ \e,n — 1) (red) sideband corresponds to the 
maximum of the narrow resonance, whereas the blue side- 
band falls into the region of the spectrum of small excita- 
tion probability, as shown in Fig. lb. The conditions for 
enhancing the red-sideband absorption while eliminating 
the carrier is therefore: 



(5) 



In summary, due to the FIT condition no absorption 
happens unless the ion is moving, and by adjusting the 
coupling Rabi frequency such that the bright resonance 
matches the red sideband, absorption accompanied by 
the loss of one phonon is made much more probable than 
absorption accompanied by gaining one phonon. The 
fact that the diffusion normally caused by \n) — s- \n) 
absorption is eliminated is most important for three- 
dimensional cooling, since the cooling limit (n) becomes 
independent of the projection of the vibrational mode on 
the laser direction 1 15[ . 



FIG. 3. Onset: (n) as a function of time calculated with 
full Monte Carlo simulation (solid line) and rate equation 
(dashed line). Parameters are Q.r = "f, Qg — 7/20, v — 7/IO, 
77 = 0.145, Ag = Ar = 2.57. Inset: Steady state distribution 

Pn- 



D. Example calculations 



EIT-cooling shows to be a robust technique even if the 
system involved deviates from the ideal 3-level scheme. 
For example, in the experiments described below, due to 
non-ideal polarizations of the laser beams, photons were 
scattered from a fourth level which introduced extra heat- 
ing. The theoretical cooling limit for this particular case 
increases by about a factor of 3, and the experiment still 
yielded cooling close the ground state |lT 



[II 



C. Cooling dynamics 



see Sec. 

Another interesting situation is when the excited state 
decays through an additional channel into a fourth level 
from which it has to be repumped with an additional 
laser. We studied the example of a single trapped IIg+ 
ion where the states \g), \r), |e) for EIT-cooling are taken 
tobetheSi/2(-F = l,7np = 1,0) andPi/2(-F' = l,mF = 
1) levels, respectively. The Si/2(F = 0, m^ = 0) level 
into which |e) decays with 1/3 probability is pumped out 
resonantly back to |e). Our model calculation is shown 
in Fig. |j. A cooling limit of (n) ~ 1 is achieved which 
is significantly below the Doppler cooling limit of about 
20 and which could help to improve the IIg+ single ion 
frequency standard [[l8[. 



A calculation starting from the full quantum mechan- 
ical master equation |16[ shows that in the Lamb-Dicke 
regime and below saturation of the cooling transition, the 
EIT-cooling process can be approximated by the same 
rate equation discussed above, only that now the appro- 
priate coefficients A± contain the quantum interference 
around A^r = A^- . The explicit form of the coefficients is 
given in p5[ . Fig. shows the cooling dynamics calcu- 
lated from a Monte-Carlo simulation of the master equa- 
tion and from the rate equation approximation. 
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FIG. 4. Onset: (n) as a function of time calculated with 
full Monte Carlo simulation for a single trapped Hg^ ion. 
Parameters are 7 = 64 MHz, !/ = 1.5 MHz, fi,. = 21 MHz, 
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0. Inset: Steady state distribution p„. 



Finally, it would certainly be also interesting to ap- 
ply EIT-cooling in cold-atom experiments, e.g. to Rb 
or Cs. In fact, the outer 77i-states of any J ^ J tran- 
sition offer a suitable level systems for its implementa- 
tion. Specifically, we give an estimate for the application 
of EIT-cooling to ®^Rb atoms in a C02-laser standing 
wave dipole trap |l3|. One main difference to an ion 
trap are that the trap frequency is generally lower and 
varies along the laser beam waist. We assumed a distri- 
bution of trap frequencies between 25 and 50 kHz, cor- 
responding to the intensity variation over one Rayleigh 
length. We applied the rate equation model to the 
(F = l,mF = 0,1) - {F' ^ l.mp = 1) 3-level sys- 
tem and included resonant repumping from F ~ 2. The 
detailed parameters and the calculated results are dis- 
played in Fig. 0. According to that estimate it should 
be possible to cool to (n) < 1, much below the Dopplcr 
cooling limit, on a time scale of ms. 



A. Levels and transitions in Ca 

The theoretical background of the method described 
above has to be modified only slightly to apply it to our 
experiment with a single trapped Ca+ ion. We imple- 
mented the EIT-scheme on the S1/2 -^ P1/2 transition. 



whose Zeeman sublevels m ~ ±1/2 constitute a four-level 
system. We denote these levels by |S,±) and |P,±), see 
Fig. |. Three of the levels, |S,±> and |P,4-), together 
with the a~^- and 7r-polarized laser beams, form an ef- 
fective three-level system of the kind considered above 
|15[ . The modifications due to the fourth level |P, — ) are 
discussed in Sec. 
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FIG. 5. Estimate for EIT-cooling of *^Rb atoms in a 
C02-laser standing wave trap: cooling limit (left) and cool- 
ing time (right) for trap frequencies of 25 .. 50 kHz. Pa- 
rameters are 7 = 6 MHz, Q^ = 1.2 MHz, Q.g = 0.1 MHz, 

ZAt- — ^9 — ^^ iVinZ, ^ Irepump — r iVlllZ, ^repump — U. 

It should be remarked that in this case the motional 
state before cooling may be outside the Lamb-Dicke 



regime, such that initially the considerations of Sec. II A 



do not apply. Instead, higher order sidebands, i.e. tran- 
sitions with |n — n'l > 1 have to be taken into account. 
In the example of Fig. |5| one finds that the absorption 
probability for cooling transitions is larger than that for 
heating transitions for many higher-order sidebands (over 
several 100 kHz), therefore starting outside the Lamb- 
Dicke regime does not seem to restrict the application of 
EIT-cooling. 



III. EIT-COOLING OF A SINGLE CA^ 
THE GROUND STATE 



ION TO 



FIG. 6. Levels and transitions in Ca used in the ex- 
periment (left side). The S1/2 to P1/2 transition is used for 
Doppler cooling and for EIT-cooling, and the scattered pho- 
tons are observed to detect the ion's quantum state. The 
narrow S1/2 to D5/2 transition serves to investigate the vi- 



brational state (see section [II C ) . Shown on the right are 
the Zeeman sublevels of the S1/2 and P1/2 states (denoted 
by |S, ±) and |P,±)) and the laser frequencies relevant for 
cooling. The solid lines labeled by tv and a"'' are used for 
EIT-cooling, with A,r = Ao-. Dashed lines indicate further 
transitions which involve the |P, — ) level. 



B. Experimental setup 

The Ca+ system seems ideal for the demonstration of 
the EIT-cooling method since in previous experiments we 
have found very small heating rates [gO| and therefore the 
cooling results can be measured with high precision. We 
will briefly introduce the experimental setup, for details 
we refer to [ po|j2l| ]. First we describe all light sources 
which are necessary for the experiment: 

For the Si /2 ^Pi/2 dipole transition at 397 nm, light 
from a Ti:Sapphire laser near 793 nm is frequency dou- 
bled in a LBO crystal inside a build-up cavity. Up 
to 5 mW of light near 397 nm are sent to the trap 
setup through a single mode UV fiber. As displayed in 
Fig. n, we use the beam deflected into -|-P* order from 
an acousto-optical modulator (AO 1) driven at 80 MHz 
for Doppler cooling. Typically, this radiation is red- 
detuned from the resonance by 20 MHz. The two beams 
for EIT-cooling, labeled a^ and tt, are derived from the 



UV fiber output as the +2"'' deflection orders of two fur- 
ther acousto-opical modulators AO 2 and AO 3 driven 
at 86 MHz and 92 MHz, respectively. We chose a de- 
tuning of Act = A^ ~ 75 MHz (natural lincwidth of the 
transition is 20 MHz). The difference in the AO 2 and 
3 drive frequencies compensates for the 12 MHz Zeeman 
splitting of the |S, ±) states in the quantization B-field of 
4.4 Gauss. The three output beams are focused onto the 
single ion in the Paul trap from different directions corre- 
sponding to their respective polarisations and the orien- 
tation of the B-field. For a reproducible and fine tuning 
of the light intensity (typically few 10 /uW in a waist of 
~ 60 fiTa) the rf-drive power of all AO's is computer- 
controlled. The rf-power, and therefore the light beams, 
can also be switched on and off within a few ^s. 
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FIG. 7. Optical setup for the light beams at 397 nm. Inset; 
Scheme of the frequencies, as used for the Doppler- and for the 
EIT-cooling. Both EIT-beams are blue detuned. The Doppler 
cooling beam is red from Ca"*" ion resonances | — , 5") -^ | — , ^) 

and l-^,^) ^|+,P). 

Light from a second TiiSapphire laser near 729 nm is 
used for excitation on the quadrupole transition from the 
Si/2 ground state to the metastable D5/2 state (sponta- 
neous lifetime 1 s). The laser frequency is stabilized to 
a bandwidth of Sv < 100 Hz using a stable high-finesse 
optical reference cavity and a Pound-Drcver-Hall servo 
circuit. The laser frequency can be tuned by an acousto- 
optical modulator, its intensity is controlled through the 
rf-drive power. With this laser we can selectively ex- 
cite vibrational sidebands of the S1/2— D5/2 transition at 
a few MHz detuning from the carrier. The relative ex- 
citation probability on the first red and the first blue 
sideband reveals information about the ion's vibrational 



quantum state as described in sect. IIIC. 

The light of grating stabilized laser diodes near 866 nm 
and 854 nm serves to pump the ion out of both 
metastable D levels. While the light of the laser at 
866 nm is turned on continuously troughout the experi- 
ment, the laser at 854 nm is also switched on and off by 
a computer-controlled acousto-optical modulator. 

Our ion trap is a 3D-quadrupole Paul trap: A ring 
with inner diameter of 1.4 mm is formed from a Mb- 
wire of 0.2 mm diameter, and two endcaps are made out 
of sharpened tips of the same material at 1.2 mm dis- 
tance pOi. The ring electrode is driven by ~ 650 V ra- 



dio frequency at 20.9 MHz, while the tips are held near 
ground potential. For the experiments described below, 
a single '*°Ca+ ion is generated by electron bombard- 
ment of a weak thermal Calcium beam. Its oscillation 
frequencies in the trap potential {vx, t^y, Vz) are (1.69, 
1.62, 3.32) MHz. Additional compensation electrodes, 
together with a differential voltage applied to the tips, 
can be used to compensate for stray static electric fields 
and to shift the ion into the center of the rf-potential. 



C. Measuring the vibrational quantum state 

The average motional quantum number rh can be mea- 
sured with high precision using laser excitation of vi- 
brational sidebands, which are spectrally well resolved 
since the linewidth of the transition, here S1/2— D5/2, 
is well below the vibrational frequencies of the trapped 
ion. The situation is schematically depicted in Fig. ra, 
where the level scheme obtained for a two-level atom in 
a harmonic trapping potential is shown together with the 
carrier (|n) — \n)) and the first sideband (|n) — \n± 1)) 
transitions. 
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FIG. 8. Principle of temperature measurement. Transi- 
tions from |S,n) to |D,n±l) are selectively laser-excited and 
the relative excitation probability is measured, vo is the bare 
atomic transition frequency. 

The excitation probabilities Pred and Pbiue on the 
\S,n) to \D,n — 1) transition (red sideband) and the 
\S,n) to \D^n + \) transition (blue sideband) arc pro- 
portional to n and ?7 -|- 1, respectively Pred drops to zero 
for |n = 0), see Fig. @. For a thermal phonon distribution 
with mean phonon number m, the probability for finding 
n phonons is Pnirh) = fn^ /(rfi -f 1)". Thus rh is directly 
related to the measured ratio Pred/ Pbiue = fh/ {rh -\- 1). 
It can be shown that this result holds for incoherent and 
coherent excitation if the phonon distribution pn{fn) is 
thermal. Our cooling results will be specified either as rh 
or by the zcro-phonon occupation pa for the three trap 
vibrational modes. 

One way to determine the cooling result is to first ex- 
cite a vibrational sideband at 729 nm and then probe the 
Si/2— Pi/2 transition at 397 nm (shelving method) [p|. 
If sideband excitation to the D5/2 state was successful, 
no fiuorescence is emitted on the S1/2— D5/2 transition. 



whereas in the other case we observe many scattered pho- 
tons. This procedure is repeated 100 times on both the 
red and the blue sideband, and the difference in their 
excitation probabihtics yields to and po- 

The second method we used is driving coherent Rabi 
oscillations on the blue sideband pq,p3|. We record 
the excitation probability Pbiue{t) as a function of 
729 nm excitation time t. Theoretically, Pbiue{t) — 
^Pni'fn) sva.^qD.nauy/'n + 1 i), where ^Rabi denotes the 
Si/2— D5/2 Rabi frequency and 77 the Lamb-Dicke factor 
p3[ . The analysis of a time evolution Puue (t) reveals the 
mean phonon number fh or can even be used to measure 
all coefficients of the phonon distribution pn(fh). 



D. Setting the power level for the EIT-beams 



E. EIT-cooling experimental procedure 

We investigated EIT cooling of the y and the z os- 
cillation at 1.62 and 3.32 MHz, respectively. The ex- 
periment proceeded in three steps: Doppler prccooling, 
EIT-cooling and finally determination of the vibrational 
quantum state, see Fig. y: 
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Applying the EIT-cooling principle to our system, we 
find that the intensity of the EIT-cr+ beam should be such 
that the Stark shift equals the frequency of the vibra- 
tional mode to be cooled. If the waist sizes of the beams 
arc known, their intensities can be calculated. There are 
also a few possible methods to experimentally determine 
the ac-Stark shift or the corresponding Rabi frequencies: 

a) A rough estimate may be gained from the fluorescence 
rate as a function of light power. Finding first the ap- 
proximate saturation power level, one can then deduce 
the necessary laser power for both EIT-beams. 

b) In the {Si/2,Pi/27D3/2} three-level system excited by 
the 866 nm and 397 nm lasers, one can measure the fluo- 
rescence rate if one of the laser frequencies is tuned over 
the resonance and the other is kept fix. Optical Bloch 
equations are then used to fit the excitation spectrum 
and determine the relevant Rabi frequencies |24|. 

c) The time constant of optical pumping between the 
|S,±) states may be used. To probe their population, 
one of the two states must be selectively excited to the 
D5/2 level. 

d) One can measure the broadening and ac-Stark shift of 
the narrow Si /2 to D5 /2 transition when the ion is simul- 
taneously illuminated with a pulse of 729 nm and one of 
the EIT-beams. 

We used the last method to determine the ac-Stark 
shift (5 as a function of power in the EIT-beams. From 
this, we could extrapolate to the power level where 
<5 — vtrap- In more detail, the ion was excited on the 
Si/2("i — 1/2) -^ D5/2(to, = 5/2) transition with a 3 /is 
pulse at 729 nm and simultaneously with a EIT-(t+ pulse. 
Although the method suffered from asymctric excitation 
line shapes and from rapid optical pumping out of the 
Si/2,"^ = 1/2 level, it lead to a good estimate for the 
right power level. For a final optimization of the EIT- 
CT+ beam power we used the cooling results for a specific 
vibrational mode |1^ . 



FIG. 9. Experimental pulse sequence. A 1.5 ms period of 
Doppler cooling is followed by a period of EIT-cooling (var- 
ied between and 7.9 ms). The cr'^'-beani is left on 50 /is 
longer than the 7r-beam to optically pump the ion into |S, -|-). 
The analysis of the vibrational state uses a pulse (of length 
t) at 729 nm for sideband excitation, followed by a Doppler 
pulse at 397 nm to probe the shelving. The laser at 866 nm 
continuously empties the D3/2 level. The whole sequence of 
20 ms duration is repeated 100 times for each sideband and 
the average excitation rate is calculated. To avoid 50 Hz mag- 
netic field fluctuations the sequence is triggered by the electric 
power line. 



(i) We Doppler precooled the ion on the S1/2 to P1/2 
transition at 397 nm (natural linewidth F ~ 20 MHz). A 
detuning of approximately -20 MHz with respect to the 
Si/2 — Pi/2 transition line was chosen. To avoid optical 
pumping into the D3/2 states, we used the repumping 
beam near 866 nm |21[| . Doppler cooling was applied for 
1.5 ms. The theoretical Doppler cooling limit on this 
transition of 0.5 mK corresponds to mean vibrational 
quantum numbers of rfi^ ~ 3 and to^ « ffi^ sa 6. The 
cooling limits reached in our experiment were higher, due 
to the fact that the simple assumption of a two level sys- 
tem in the determination of the Doppler limit does not 
hold in our case and because the 12 MHz Zeeman split- 
ting of the Si/2 state docs not permit optimum detun- 
ing for all transitions between the substates. We experi- 
mentally determined the mean excitation numbers after 
Doppler cooling to be m^ = 6.5(1.0) and m** = 16(2). 

(ii) After Doppler cooling, we applied the two EIT- 
beams, 397nm-cr^" and 397nm-7r for up to 7.9 ms. A 
1.8 ms pulse duration was found to be sufficient, longer 
cooling time did not lower the mean phonon number (see 
Sec. IV B where we will discuss the cooling time). The 
power levels were set as discussed above in Sec. HID. 
For our setup with laser beam waists of ~ 50 /xm, a laser 
power of < 50 fiW was used for the cr+-bcam. 



The fc- vectors of the two EIT-cooling beams enclosed an 
angle of 125° and illuminated the ion in such a way that 



their difference Afc had a component along all trap axes 
{{(j)x,4>y,4>z) = (66°, 71°, 31°), where (/)i denotes the angle 
between Afc and the respective trap axis). As a result, all 
vibrational directions were cooled. Since the vacuum re- 
cipient did not allow for two orthogonal beams, we could 
not realize the ideal constellation of Fig. g: While the (T+- 
beam was parallel to the quantization axis as desired, the 
TT-beam was not orthogonal to it and therefore had a a~- 
componcnt. The extra excitation due to this non-ideal 
constellation is indicated in Fig. || by a dashed line. It 
results in a slightly higher final mean phonon number. A 
detailed discussion is found in Ref. [^ . 

(iii) We analyzed the vibrational state after EIT- 
cooling by spectroscopy on the S1/2 — *■ D5 /2 qu adrupole 



transition at 729 nm, as described in Sec. IIIC, Finally 



a diode laser at 854 nm served to rcpump the ion from 
the D5/2 to the S1/2 level. 



IV. EIT-COOLING RESULTS 

A. Cooling results for a single mode of vibration 

With the ac-Stark shift 6 set to the frequency of the 
radial y-modc, we monitored the vibrational state after 
EIT-cooling by exciting the blue sideband of the |S, +) -^ 
D5/2('^ = +5/2) transition with a 729 nm pulse and then 
measuring the |S, +) level occupation as a function of the 
pulse length t [^,0. The observed Rabi-oscillations 
were subsequently fitted to determine the mean vibra- 
tional occupation number to^ |22|, see Fig. |o[ The low- 
est mean vibrational number m** = 0.18 observed corre- 
sponds to a 84% ground state probability. We repeated 
this experiment on the z-mode at 3.3 MHz after hav- 
ing increased the intensity of the cr^-bcam to adjust S. 
For this mode, a minimum mean vibrational number of 
fh^ =0.1 was obtained, corresponding to a 90% ground 
state probability. 
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FIG. 10. Rabi oscillations excited on the blue sideband of the radial y-mode (left) and the axial z-mode (right) of a singl e 
ion. From the theoretical curves for Puucij), mean phonon numbers of m^ = 0.18 and m^ — 0.1 are determined (see Sec. IIIC). 



We found the cooling results largely independent of 
the intensity of the 7r-beam as long as it is much smaller 
than the a~^ intensity. In our experiment the intensity 
ratio was la / I-k — 100 and we varied the intensity of the 
TT-beam by a factor of 4, with no observable effect on the 
final TO. 



B. Cooling dynamics 
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By determining the dependence of the mean vibra- 
tional quantum number on the EIT cooling time, we 
measured the cooling time constant for the y-mode to 
be 250 /iS, as shown in Fig. |l^. 



FIG. 11. Mean vibrational quantum number fa" versus 
EIT cooling time. The insets show Rabi oscillations excited on 
the upper motional sideband of the |S, +) — > 05/2(771 = +5/2) 
transition, after 0.4 ms (left) and 7.9 ms (right) of EIT cool- 
ing. A thermal distribution is fitted to the data to determine 
fh" in both cases. 



C. Cooling of two modes 

To show that the EIT method is suitable to simul- 
taneously cool several vibrational modes with signifi- 
cantly different frequencies of oscillation, we chose the 
axial z-mode at 3.3 MHz, and the radial y-niode at 
1.62 MHz. The intensity of the iT+-beam was set such 
that the ac-Stark shift was roughly halfway between the 
two mode frequencies. Again we applied the EIT cool- 
ing beams for 7.9 ms after Doppler cooling. This time 
we determined the final m^'^ by comparing the excita- 
tion probability on the red and the blue sidebands of the 
Si/2("i = 1/2) -^ D5/2(?Ti = 5/2) transition. The result 
is shown in Fig. |lj. We find both modes cooled deeply 
inside the Lamb-Dicke regime [rj^/Wi ~ 0.02 <C 1), with 
Pq = 58% and pg = 74% ground state probability. 
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FIG. 12. EIT-cooling of two modes at 1.6 MHz and 
3.2 MHz simultaneously. From the sideband excitation rate 
after cooling we deduce a ground state occupation number 
of 73% for the axial mode (3.2 MHz) and 58% for the radial 
mode (1.6 MHz). 



V. EIT-COOLING OF LINEAR ION STRINGS 

Since EIT-cooling allows simultaneous cooling of sev- 
eral modes at different frequencies, it seems to be partic- 
ularly suited for ion strings in linear traps. The frequen- 
cies of the axial vibrational eigenmodes of a linear string 
have been calculated (25U2q] and measured p^. For a 
10-ion string trapped in a linear trap with a center-of- 
mass axial frequency of 0.7 MHz, the closest inter-ion 
spacing is found to be 3.0 /xm (4.5 /.tm for N=5). The 



axial vibration frequencies are 0.7 MHz, 1.22 MHz, .. 
4.6 MHz. The radial trap frequencies in a linear trap 
must be made sufficiently high in order to prevent a 
transition from the linear configuration of N ions to a 
zig-zag configuration. It was estimated that this tran- 
sition occurs at Vradjvax ^ 0.73iV°-*^, thus the radial 
trap frequency must exceed 3.7 MHz for 10 ions (2 MHz 
for N=5) [^. Typically, the radial frequency is cho- 
sen higher. The linear ion trap experiment at Innsbruck 
uses Vrad > 4 MHz, the Be+ experiments at NIST |2|] 
have Vrad > 20 MHz. Apart from the purely axial modes 
there exist 2N radial modes. Their frequencies arrange 
in a band below Vrad which overlaps with the axial mode 
frequencies. 

We now estimate the performance of EIT-cooling for 
a 10-ion string. The result is displayed in Fig. |l3|: For a 
10-ion string indeed all 3N vibrational modes are cooled 
to a mean phonon number m below one. This is promis- 
ing for the application of cold ion strings for quantum 
information processing [§|-||jlj . As discussed in the intro- 
duction, it is required that all modes which couple to the 
laser light (spectator modes) must be cooled well into the 
Lamb-Dicke regime. The reason for that is that thermally 
excited spectator modes cause a blurring of the Rabi fre- 
quency CI which disturbs the precision of quantum oper- 
ations. For the case of 3N — 1 spectator modes i with 
to' and Lamb-Dicke factors 77^, the initial state is a mixed 
state. Each time the experiment will happen with slightly 
different initial conditions and the Rabi frequency for car- 
rier or sideband transitions will differ. The relative blur- 
ring reads like Af^/rj = ^^ rjf fh' [fh^ + l)/(3iV - 1) 
(see equ. 126 in ref. [p8|) which we find as small as 3- lO""' 
for our specific example. The maximum number of Rabi 
oscillations will then be -- (AO/0)"i -- 2000. 



cT) 


0.7 







F 


U.b 


-1 




c 


U.b 


r 







4 


c 






r 


0.3 


Q. 




r 


0.2 


(fl 




<u 


0.1 




12 3 4 

vibr. mode frequency (MHz) 

FIG. 13. EIT-cooling of the axial modes of a lin- 
ear string based on the S1/2 - P1/2 transition with 
0.^= 30 MHz, n^= 0.5 MHz, Vp= 20 MHz, and a detuning 
of Ao- = A.n-= 75 MHz. The axial trap frequency is 0.7 MHz. 
For the calculation we have chosen the light intensity such 
that the bright state is ac-Stark shifted by ~ 3 MHz (see in- 
set). The bright resonance with a width of ~ 0.5 MHz leads 
to cooling for all axial modes. The mean phonon numbers 
(black dots) of all axial modes are plotted versus the mode 
frequencies. 



VI. CONCLUSION 



In summary we have presented EIT-cooling, a novel 
ground state cooling method which uses electromagneti- 
cally induced transparency on two coupled dipole transi- 
tions. We have demonstrated its successful application to 
a single Ca+ ion in a Paul trap. We have also shown sev- 
eral estimates for its implementation in various trapped- 
ion and trapped- atom experiments, taking into account 
different parameter regimes as well as deviations from the 
ideal case such as extra levels and transitions. The EIT- 
cooling method, according to these estimates and consid- 
ering its experimental simplicity, bears the potential to 
become a standard tool for the preparation of ultracold 
atoms and ions. A particularly important application 
will be cooling of an ion string in a linear trap, as a step 
in preparing such a string for coherent manipulations in 
quantum information processing. 
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